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We present results for the spectrum of light and strange mesons on configurations with two 
flavors of mass-degenerate Chirally Improved sea quarks. The calculations are performed on seven 
ensembles of lattice size 16 3 x 32 at three different gauge couplings and with pion masses ranging 
from 250 to 600 MeV. To reliably extract excited states, we use the variational method with an 
interpolator basis containing both gaussian and derivative quark sources. Both conventional and 
exotic channels up to spin 2 are considered. Strange quarks are treated within the partially quenched 
approximation. For kaons we investigate the mixing of interpolating fields corresponding to definite 
C-parity in the SU(3) limit. This enlarged basis allows for an improved determination of the low- 
lying kaon spectrum. In addition to masses we also extract the ratio of the pseudoscalar decay 
constants of the kaon and pion and obtain Fk/F^ = 1.215(41). The results presented here include 
some ensembles from previous publications and the corresponding results supersede the previously 
published values. 

PACS numbers: 11.15.Ha, 12.38.Gc 



I. INTRODUCTION 

Considering only strong decays, with the exception of 
the pion and the proton all hadrons are resonances, em- 
bedded in a continuous spectrum. In lattice calculations 
we can only determine discrete energy levels, with spac- 
ings 0(1/L) related to the spatial extent L of the studied 
lattice volume. When disregarding the fermion vacuum 
in the so-called quenched simulations energy levels can be 
related directly to hadron excitations. In dynamical situ- 
ations the energy levels are denser close to resonances and 
they are influenced by coupled open hadronic scattering 
channels. Although in principle the Euclidean correla- 
tor of any hadron interpolator with the correct quantum 
numbers should feel these scattering channels, in actual 
calculations there is little, if any, trace of it [l|, unless 
such multi-hadron interpolators are included explicitly in 
the set of operators. However, inclusion of those is costly, 
since it involves disconnected contributions. In actual 
calculations efficient but demanding all-to-all propagator 
methods are used 0-0] • 

In recent years much effort has been invested into de- 
veloping methods for determining the lowest energy levels 
for hadron correlators. In fil. I2L Tsl— T]~3l| meson excitations 
have been studied in a dynamical quark background with 
a variety of quarks species, interpolators and extraction 
methods. A central technique employed was the varia- 
tional method [3, [l5| where one finds the energy levels 
by diagonalization of cross-correlations of a (hopefully) 
sufficiently large set of interpolators which allows for a 
good overlap with the relevant hadron states. 

In continuum quantum field theory there has been 
recent progress in investigations of mesons using 
Schwingcr-Dyson equations and the Bethc-Salpeter equa- 



tion as well as effective field theories (see for examples 

Refs. mil). 

Starting with [23| we have been determining hadron 
ground states and low excited states in a framework 
of simulations with two light dynamical quarks. The 
fcrmionic action used is the so-called Chirally Improved 
(CI) action [24|, [25[, an approximate solution to the 
Ginsparg- Wilson relation for fermions obeying chiral 
symmetry in a lattice form. The strange quarks have 
been incorporated in the valence sector only. In [l[ re- 
sults based on three ensembles at three different gauge 
couplings but with only one quark mass for each coupling 
have been presented. We have meanwhile significantly 
extended the statistics and also the number of ensem- 
bles. Here, we present our results for the meson sector 
based on the final set of seven ensembles at three gauge 
couplings and two or three quark mass values at each. 
This allows an extrapolation towards the physical point. 
Previously published results are generally confirmed, al- 
though in some cases we observe new behavior related to 
new symmetry considerations. Some results have been 
presented already in [2?jj . 

Following the presentation of the action and the pa- 
rameters of the gauge configuration ensembles in Sec. |TT] 
we discuss scale setting, decay constants and the quark 
mass in Sec. HID The interpolators used for the meson 
fields in the variational analysis are discussed in Sec. IIVI 
and tabulated in the appendix. The main parts are Sees 
IVl and fVTl where results for the mesons are presented. 
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II. ACTION AND SIMULATION 

A. Fermion action and gauge action 

In our study the fermions are represented by the Chi- 
rally Improved Dirac operator Dei @, HH • This is an 
approximate solution of the Ginsparg- Wilson equation 
and results from a general ansatz for the Dirac operator, 
namely an expansion of the form 

16 

D = m l + D C i , D cl (n, m) = £ c^, (17) T t , (1) 

i=l 

where the sum runs over all 16 elements I\ of the Clifford 
algebra and the coefficients elm were fit by minimizing 
the violation of the Ginsparg- Wilson equation. It in- 
cludes paths up to a maximum length of 4 lattice units. 
The paths and coefficients used are found in the appendix 
of j23[. We used the same 19 coefficients for all ensem- 
bles, modifying only the diagonal mass term in order to 
account for the additive mass renormalization. For that 
reason the values of the bare mass parameter mo given in 
Table U are negative. Thus the actual (unrenormalized) 
mass is given by the values ttiawi determined from the 
axial Ward identity. 

For further improvement of the fermion action one level 
of stout smearing of the gauge fields [13 was included in 
its definition. The parameters arc adjusted such that the 
value of the plaquette is maximized (p = 0.165 following 
[Ij). For the pure gauge field part of the action we use 
the tadpole-improved Luscher-Weisz gauge action [28[. 
For a given gauge coupling we used the same assumed 
plaquette value for the different values of the bare quark 
mass parameter. 

B. Lattice ensembles 

The analysis presented here is based on seven ensem- 
bles of configurations for lattice size 16 3 x 32. These sub- 
stantially extend (by a factor of three) the data base of 
[l|, [23| . A summary of the notation and some parameters 
of these ensembles is given in Table fl] 

The notation for the couplings follows (23|, where all 
parameters of the fermion action are detailed. For each 
value of the gauge coupling we have two or three values 
of the quark mass parameter. Following equilibration ev- 
ery 5th configuration has been selected for analysis. Fur- 
ther details on the updating HMC-method and statistical 
checks for equilibration have been discussed in [23| . 

From the values of m x L we expect non-negligible finite 
size effects for the three ensembles with smallest quark 
mass, A66, B70 and C77. Discretization effects have been 
discussed in the quenched simulations, where for the used 
action only small 0(a 2 ) corrections have been identified 
[29l ]. In order to confirm this for the dynamical simula- 
tion we would have to perform our study at several lattice 
spacings and volumes, which is not possible based on the 



set 


Plw 


m 


m 3 


configs 




A50 


4.70 


-0.050 


-0.020 


200 


6.4 


A66 


4.70 


-0.066 


-0.012 


200 


2.7 


B60 


4.65 


-0.060 


-0.015 


300 


5.7 


B70 


4.65 


-0.070 


-0.011 


200 


3.4 


C64 


4.58 


-0.064 


-0.020 


200 


6.7 


C72 


4.58 


-0.072 


-0.019 


200 


5.1 


C77 


4.58 


-0.077 


-0.022 


300 


3.7 



TABLE I. Parameters of the simulation: We used several 
ensembles with different gauge couplings /3lw and/or light 
quark mass parameters m = 0. We also show the strange 
quark mass parameter m s , the number of configurations ana- 
lyzed and the physical extent of the spatial volume multiplied 
with the pion mass. 



given ensembles and statistics. Studies with a larger vol- 
ume (24 3 x 48) with linear size 0(3.6 fm) are in progress. 



III. SCALE AND LOW ENERGY PARAMETERS 
A. Scale 

In our earlier work [l], [23[ we had analyzed configu- 
rations at one quark mass parameter for three values of 
the gauge coupling. There, we used the lattice spacing 
derived from the static potential with a Sommer param- 
eter ro = 0.48 fm. Now we have two or three quark mass 
parameters for each gauge coupling and can attempt an 
extrapolation to the physical point or the chiral limit. 
The latter extrapolation would be relevant for the pa- 
rameters of Chiral Perturbation Theory (ChPT), which 
we will not attempt to extract here. 

We use two approaches to set the scale. In the first 
one we determine y = a/rg from the static potential sep- 
arately for each ensemble, as discussed in [2j]. We then 
study the dependence of this quantity on the measured 
values of x = (am^) 2 (cf., Fig. [1]). The physical values 
are obtained along 



For each of the three gauge couplings we then perform 
a linear fit and obtain the physical value where the ex- 
trapolations intersect Eq. (|2|) with m^r® = 137 MeV x 
0.48 fm = 0.3332. (Wc use the average of charged and 
neutral pion masses.) From this one reads off the lat- 
tice spacing a. Table [TT| gives the resulting value in the 
row labeled (TT-,ro) p hys- The value in the chiral limit is 
obtained as usual from a/ro where am w = 0. 

The other approach is to replace y = a/ro by mass 
values like ampf or am p . Since the p is unstable for small 
enough pion mass, there will be threshold effects. In our 
parameter range we find no coupling to the (p-wave) irir 
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FIG. 1. (color online). Setting the scale with the Sommer 
parameter and the pion mass as input at the physical point. 
The green (long-dashed) line is the curve Eq. ([2]). The solid 
and short-dashed lines represent the extrapolation of our lat- 
tice data. Their intersections with the green line define the 
lattice constants a. 





A 


B 


C 


(tt, ra) P hys 


0.1324(11) 


0.1366(15) 


0.1398(14) 


(l", ro)chiral 


0.1314(12) 


0.1356(17) 


0.1387(15) 


(■K,p)ph V s 


0.1330(44) 


0.1378(50) 


0.1400(29) 



TABLE II. Lattice spacing in physical units derived for en- 
sembles of type A, B, C (cf., Table[T| by the methods discussed 
in the text. 



sector yet and a linear extrapolation intersecting with 
y = ^/xirip/m^ gives the values of the lattice spacing in 
Table [TT1 compatible with the results of the first method, 
but with larger errors. 

Throughout this presentation we will use the values 
obtained from the definition denoted by (w, ro) p h ys in 
Table HU 



B. Setting the strange quark mass 

In this two-flavor simulation we use the partial quench- 
ing approximation to access the strange hadron spec- 
trum, i.e., we consider the strange quark as a valence 
quark only. In view of results with full strange quark 
dynamics (e.g., [30| ) we find, at least for the ground 
states, no noticeable difference in the mass range con- 
sidered here. In each ensemble the strange quark mass 
parameter m s is set by identifying our result for the f2 
baryon positive parity ground state energy level with the 
physical 17(1672). These parameters are found in Table 

in ' 

For this definition we use rQ^ exp = 0.48 fm in each en- 
semble, differing from the (in Subsec. IIII A[) discussed 



Set 


a 


a m-K 




am AWI 


m AWI 




[fm] 




[MeV] 




[MeV] 


A50 


0.1324(11) 


0.3997(14) 


596(5) 


0.03027(8) 


45(1) 


A66 


0.1324(11) 


0.1710(48) 


255(7) 


0.00589(40) 


9(1) 


B60 


0.1366(15) 


0.3568(15) 


516(6) 


0.02356(13) 


34(1) 


B70 


0.1366(15) 


0.2111(38) 


305(6) 


0.00836(23) 


12(1) 


C64 


0.1398(14) 


0.4163(18) 


588(6) 


0.02995(20) 


42(1) 


C72 


0.1398(14) 


0.3196(18) 


451(5) 


0.01728(16) 


24(1) 


C77 


0.1398(14) 


0.2340(27) 


330(5) 


0.01054(19) 


15(1) 



TABLE III. Pion masses and quark AWI-masses for the dif- 
ferent sets of gauge configurations. 



method to set the overall scale. Since the two different 
definitions agree at physical pion masses, this method is 
consistent at the physical point, but results have to be 
taken with care at unphysically large pion masses. 



C. AWI mass 

The so-called axial Ward identity (AWI) mass (or 
PC AC mass) is determined from the asymptotic (i.e., 
plateau of the) ratio of the unrenormalized correlators 



2 mAwi = 



c A (0|cVl+(p = 0,t) X(0)\0) 
c P (0|P-(p = 0,t)X(0)|0) 



(3) 



where P~ = d^u, = CZ7475U, and X is an interpola- 
tor with the quantum numbers of the pion, usually P + 
or A + . The constants ca{s) and cp(s) denote the lat- 
tice factors relating the smeared interpolators to the lat- 
tice pointlike interpolators (not to be confused with the 
renormalization constants Z relating lattice point oper- 
ators to the continuum renormalization scheme). They 
are obtained from the ratio of correlators from smeared 
to point sources [23^ . 

The relation to the renormalized quark mass needs the 
renormalization factors for the pseudoscalar and axial 
currents, 



m 



(r) 



Z A 



■ m awi 



(4) 



Table Hill gives the values of wawi and for the ensem- 
bles studies. (Values for the renormalization constants 
have been derived in |3lL l32T|.) 



D. Decay constants 

The pseudoscalar decay constant describes the cou- 
pling to weak decays. It can be extracted from the 
asymptotic behavior of the correlation between the pseu- 
doscalar or the time components of the axial interpola- 
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0.20 



FIG. 2. The ratio of Fk/F-x is plotted against m\ (in dimen- 
sionless units) for each set of gauge configurations. The full 
black line is a fit of the data using the relevant expressions for 
numerator and denominator; the shaded area indicates the er- 
ror band. The magenta cross indicates the experimental value 

m. 



tors. 



c\ Z 2 A (A-(j> = 0, t) A+(0)) ~ m n Fl e 
The coefficient then gives 



Fjr 



2m A wi cp Zaa —t , 



(5) 



(6) 



and equivalcntly for the kaon Fk ■ 

The dependence of the pion decay constant on the 
quark mass can be described by chiral perturbation the- 
ory. Up to 1-loop order one finds [33j 



F-x — F\rJ 



2Sr 



16^^,0 



In 



2 So 



(7) 



Here, F^fi and So refer to the pion decay constant and 
the quark condensate in the chiral limit m — > and A4 
is a low energy constant. The corresponding expressions 
including the 2-loop order can be found in [3J, l35| . 

The renormalization factor Za cancels in the ratio 
Fk/F-k- We show this ratio in Fig. [2] where we assume 
a lattice spacing of 0.135 fm (the average of our values 
for the scheme (jr,ro)phys) and a physical pion mass of 
139.57 McV. The extrapolation of our data to that point 
gives 

F K /F W = 1.215(41) . (8) 
which fully covers the experimental value 1. 197(9) (36j. 



IV. ANALYSIS METHOD AND MESON 
INTERPOLATORS 

Given interpolating operators Om with the quantum 
numbers of a hadron, the correlation function of such op- 
erators separated by some Euclidean time distance pro- 
vides the energy spectrum, 



(O M (t)OU0)) = y2(0 M \n){n\Ol)e 



-E n t 



(9) 



The asymptotic exponential decay, however, gives just 
the ground state energy in that channel. On finite lat- 
tices, depending on parameters like size and lattice spac- 
ing, this may be related either to a single meson or to 
meson scattering states. For the study of scattering and 
of higher lying mesons it is imperative to find also the 
excited energy levels. 

An efficient tool for this is the so-called variational 
analysis [l4|, EH H3j- Using several interpolators with 
the correct quantum numbers, one diagonalizcs the cross- 
correlation matrix of these, using the generalized eigen- 
value formulation 



C!y(t) 

C(t)v k (t,t ) 



(Oi(i)Oj(O)) , 
Xk(t,t )C(t )v k (t,to) 



(10) 



If the set of interpolators is large enough, then one ex- 
pects that the eigenvectors approach the eigenstates of 
the system. In fact, the eigenvectors act as a fingerprint 
of the states and should remain stable over the considered 
window of t- values. In such a window the eigenvalues 
decay exponentially, approximating the desired eigenen- 
crgics, 

A fc (Mo)oce- (t -* o)Bfc (l + O{e-^~ t0 ^ AE ^ . (11) 

Here, depending on t and to the value of AE k denotes 
either the difference to the first neglected energy level 
(for t < t < 2t ) or to the nearest energy level (for a 
careful discussion see It was also demonstrated, 

that even ghost states can be identified with this type of 
analysis [38l | . 

A possible systematic influence comes from choosing 
to in the variational method and the fit range for the 
generalized eigenvalues. We use to = 1 throughout. In 
principle, the impact of that choice can be estimated by 
choosing several values of to and varying the fit range. 
For the final fit one should then choose a window where 
this impact is negligible. However, in practice the cor- 
responding choices are restricted by the given signal-to- 
noise ratio for coarse lattices and weak signals. In the 
actual analysis one determines the window from a combi- 
nation of indicators, ranging from effective energy values 
to approximate constancy of the corresponding eigenvec- 
tors. The energy levels then result from an exponential 
fit to the eigenvalues over that window. In some cases 
a second exponential is used in these fits to allow for a 
small admixture of higher energy states. 
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Various techniques have been suggested to construct 
interpolators. In [39j we introduced lattice operators 
based on smeared quarks. Combining differently smeared 
quarks, also including covariant derivatives [40l - l42l ]. sev- 
eral meson and baryon energy levels could be determined 
in the quenched [43|, |44[ and dynamical case [l[ ■ 

The interpolators are constructed by hypercubic 
(HYP)-smearing (45l - l47j the time slice gauge variables, 
i.e., smearing only the spatial links in each time slice 
0. Based on these gauge variables the quark sources are 
smeared with the covariant Jabobi smearing [48|, |49[ 



K 



S K ,K = "So , 

n=0 
3 

H(n, m) = Y (Uj (ft, 0)s(n + 



(12) 
(13) 



3=1 



+ Uj ( ft - j , 5 ( n - j, m) , 



where So denotes the point source. The parameters K 
and k are adjusted to obtain gaussian-like shapes of the 
sources [44[ with different smearing widths. In the defi- 
nitions of the operators we denote the smearing types by 
n and w (narrow and wide) and by dk for the derivative 
in spatial direction k. The widths of the sources do not 
exactly agree for the various ensembles (which would be 
dependent on the definition of the scale anyway.) How- 
ever, the width of the narrow source is in the range 0.2 
to 0.3 fm and the width of the wide source is in the range 
0.4 to 0.6 fm. 

The derivative sources Sg k have been constructed nu- 
merically by applying the covariant difference operators 
on the wide source, S w , see [42|, This corresponds to 
an asymmetric definition of the interpolators. If S^S 2 
denote gaussian smearing operators and t5 the deriva- 
tive acting to the right, then our operators (involving 
one derivative) have the structure 



instead of 



o = i>(s 1 rs 2 ~3 ± < ds 2 ts 1 )^ (u) 



o = ^(s 1 rl$s 2 ±s 2< Drs 1 )ij , (15) 



where the "±" symmetrization ensures a good C-parity 
quantum number. Following Eq. (|15p . some interpolators 
(with S 1 = S 2 ) are identical to zero after partial integra- 
tion. The operator Eq. (|14p is in general non-vanishing 
even if S 1 = S 2 , since [D,S] 7^ 0. This commutator 
can be seen as introducing additional pieces of paths in 
the combined smearing operator, which means changed 



1 Notice that the Dirac operator already contains one level of stout 
smearing. We use these stout smeared gauge links and apply 




0.2 
m 2 [GeV 2 ] 



FIG. 3. _+ (top): only the first excitation is shown, the 
ground state pion mass squared defines the abscissa. ++ 
(bottom): observed ground state and first excitation. 



additional smearing to construct the sources. 



weights of the existing paths and a few new paths. Nu- 
merically, we find that the corresponding correlators are 
of the same magnitude as others and yield consistent sig- 
nals. Hence, this asymmetric definition enlarges effec- 
tively the basis of operators to some extent. In particu- 
lar some exotic channels can be accessed this way already 
with fewer derivatives. 

In Appendix \X\ we list all meson interpolators used in 
our study, ordered according to their spin and parity. The 
tables differ from those in [l| since we here account for the 
approximate symmetry under C-parity of strange mesons 
and construct the interpolators accordingly. Monitoring 
the eigenvectors in the variational method allows for in- 
sights in approximate C-parities of various strange me- 
son states, and furthermore in the breaking of C-parity 
of strange mesons when approaching the physical pion 
mass. 
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V. ISOVECTOR LIGHT MESONS 

The energy levels are obtained from exponential fits to 
the eigenvalues in a range of t- values where the eigenval- 
ues and eigenvectors are compatible with plateau behav- 
ior. Typically that plateau extends from t = 2 or 3 up 
to t = 6 to 12. In some cases the eigenvalues are close to 
each other and their order changes from one timeslice to 
another and also changes randomly over the set of con- 
figurations. This complicates the exponential fits to the 
eigenvalues and the automatic attribution of the eigen- 
vectors to physical eigenstates. In such situations we 
use scalar products of eigenvectors at a given timeslice 
with the eigenvectors at the preceding timeslice to sort 
the eigenvalues according to their corresponding physical 
states. This procedure becomes more reliable towards 
finer lattice spacings. For subsets of configurations (in 
the jackknife analysis) the eigenvectors are contracted 
with the average of the vectors at the same timeslice. 

All masses are extrapolated towards the physical point 
as a function of the pion ground state mass squared. In 
the plots we also show the corresponding one a error band 
(dashed curves). The number of energy levels shown is 
always less than the number of interpolators chosen for 
the diagonalization. The x 2 P er degree of freedom for the 
chiral fits of all energy levels are collected in the Tables 

ED and lxrmixrvl 




FIG. 4. Eigenvector components of the ground state of the 
light scalar meson channel (ao) of ensembles A50, C77 and 
A66 (top to bottom). Interpolator (4) (only gaussian sources) 
dominates, while contributions of the other interpolators (one 
or two derivatives) is found to be particularly relevant at 
heavy pions. Nevertheless, the eigenvectors are very similar 
over the whole range of pion masses (600 to 250 MeV) and 
only evolve smoothly. 



Scalars 



l ~(7r): For the first excitation in the pion channel 
(see Fig. [3]), the set of operators (1,2,17) is used in all 
ensembles. The corresponding effective mass plateaus 
arc rather short, increasing the uncertainty of the ex- 
tracted mass. Due to the finitcness of the lattice, the 
back-running pion limits the possible fit range for the 
first excitation [2^, |42|, [5l| , in particular at small pion 
masses. Nevertheless, masses can be extracted and the 
chiral extrapolation hits the experimental 7r(1300) within 
la. 



2 
1.5 

6 i 

0.5 



o Ground state 
□ First excitation 



a (A66) 



FIG. 5. Effective mass plateaus for the light scalar meson 
channel (ao) of ensemble A66, ground state and first excita- 
tion. 



0++(a ): In Q three (A50, B70 and C77) of the seven 
ensembles have been analyzed, with less statistics than 
in the present work. Partially quenched data was used 
to argue that the signal in the ++ channel probably has 
significant contributions from the S-wave scattering state 
7r?72 • In the present work we analyze only fully dynamical 
data (except for the strange sector). Our results are now 
compatible with the experimental ground state a(980) 
within 1 <T and with the first excitation a(1450) within 
2 a (sec Fig. [3]). However, the channel still poses some 
difficulties. The plateau is rather short and there remains 
some ambiguity in choosing the fit range, leading to a 
systematic error. In addition, the results depend on the 
chosen set of interpolators. We show results from subsets 
of (1,4,10,12,13). In ensemble B60, the excitation signal 



was not good enough to be fitted. The extrapolations of 
the ground state levels agree for the different choices of 
interpolators. 

However, in particular the ground state energy level of 
ensemble A66 deviates when changing the set of inter- 
polators. The result becomes unexpectedly light, most 
pronounced in the case of the set (10,12,13), though the 
corresponding effective mass plateaus look stable. In- 
deed, this point lies below the (theoretical) wrfe thresh- 
old and could indicate a scattering state signal. It also 
could signal a severe finite size effect for this case in A66; 
this could be clarified only by increasing the lattice vol- 
ume. Nevertheless, except for this point, the results are 
compatible with the experimental states. 

In Fig. 0] we show the eigenvectors for the ground state 
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for three ensembles covering the whole range of pion 
masses presented. They are quite consistent with each 
other and not supporting the notion of a change in the 
physics of the ground state over that range. Fig. [5] shows 
the effective masses of ground state and first excited en- 
ergy level for the ensemble with smallest pion mass (A66). 

There are studies for the finite size dependence of the 
lowest energy level in this channel based on unitarizcd 
chiral perturbation theory However, at the mo- 

ment our values are not precise enough to decide on these 
grounds on properties of the ag. Also it may be necessary 
to include meson-meson interpolators in a more detailed 
study. Simulations to address finite-size effects are cur- 
rently in progress and the discussion of this ongoing effort 
is beyond the scope of the current publication. 



B. Vectors 



1 ++ (ai ) : The signal in the pseudovector meson channels 
is usually bad compared to the pion and the p channels. 
Nevertheless, the ground state and a first excitation can 
be identified. The ground state is extracted using the sin- 
gle interpolator (1). For the first excitation the set has 
to be optimized in each ensemble separately, taking sub- 
sets of three interpolators out of (1,2,4,13,15,17). Some 
of the plateaus tend to shift towards smaller masses at 
large time separations. However, as far as possible, long 
fit ranges are chosen. The chiral extrapolations hit the 
experimental ai(1260) and the ai(1640) within error bars 
(see Fig. 

I" 1 (t>i): In the 1 H channel, the ground state plateau 
is more stable than in its positive C-parity partner chan- 
nel (ai). Using the single interpolator (6), a mass with 
comparatively small error bar is obtained. The chiral ex- 
trapolation comes out too high, missing the experimental 
&i(1235) by more than 2a (see Fig. [6]). 



l~~(p): The /o(770) comes out nicely as usual (see 
Fig. [6]). The first and second excitation are extracted 
using the set (1,8,12,17,22), where the second excita- 
tion is not stable in A66. These excitations are very 
close to one another, making the chiral extrapolations 
less reliable. The pattern of energy levels would allow a 
crossover of eigenstates but the eigenvectors do not con- 
firm this. Therefore, we extrapolate the results to the 
physical point according to the naively assumed level or- 
dering, neglecting a possible crossover. The results are 
compatible with the experimental p(1450) and p(1570 or 
1700) within error bars (for a discussion on the latter 
excitation see (36^. 

We find no obvious indication for a coupled mr P-wave 
channel. As discussed earlier 0, [52[ this may be due to 
weak coupling. By including two pion interpolators one 
can derive a scattering phase shift from the modification 
of the observed energy levels close to the resonance (see, 
e.g. Hl|). Such a study needs inclusions of disconnected 
graphs, which are not accessible to us: The necessary 
propagator calculation is numerically too costly for CI 
fcrmions. 

1 (7Ti): The quantum numbers 1 cannot be ob- 
tained with isotropic quark sources only. Thus, this 
channel is not accessible by simple quark models, and 
it is commonly referred to as exotic. Due to the weak 
signal, the set of operators has to be optimized in each 
ensemble separately, taking one or two interpolators of 
(9,11,14,16,21,24). This way a mass value can be ex- 
tracted only with comparatively large statistical uncer- 
tainty. The chiral extrapolation hits the experimental 
7Ti(1400), but is also compatible with the 7i"i(1600) (see 
Fig. [5]). In some of the ensembles we get the best signal 
using interpolators which are nonzero only due to the 
definition in Eq. [H] and discussed there. This may be 
related to the "exotic" property of this channel. 



C. Tensors 

The continuum representation for spin 2 contributes to 
the irreducible representations and E on the lattice. 
These interpolators are orthogonal, thus masses can be 
extracted in each of them separately. In the continuum 
limit, the results should agree, however, at finite lattice 
spacings they can show different discretization effects. 
We extract the energy levels separately and compare the 
corresponding chiral extrapolations. 

2 (P2): In many of the spin 2 channels the signal 
is weak and fits can be performed only for some of the 
seven ensembles. In particular this is the case in the 2 
channel (see Fig. top and middle). We use the single 
interpolator (2) in T2 and also (2) in E. The effective 
masses are noisy, the fitted plateaus are rather short, 
with only 2 d.o.f. in the fits. Nevertheless, the chiral 
extrapolations of the T2 and E ground state masses agree 
with each other and also with the experimental ,02(1940) 
mass. Hence, our results are compatible with this state, 
which is omitted from the summary table of (36j . 

2 '"(^2): In the 2 h channel (Fig. [JJ bottom), interpo- 
lator (6) is applied in T2. The extrapolation to the phys- 
ical point is compatible with the experimental 7^(1670) 
(within 1 resp. 1.5a). The signal for representation E 
(not shown) is too weak to be reliable. 

2 H : We studied this channel for completeness but the 
signals were inconclusive and did not allow to extract an 
energy level. 

2 ++ (a2): In the 2 ++ channel (Fig. [5]), we use interpo- 
lator (2) in T 2 and (2) (respectively (6) for A66) in E. 
Some of the plateaus are unexpectedly light, however, 
that might be statistical fluctuation. The chiral extrap- 
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olations of the T2 and E ground state masses agree and 
both match the experimental 02 (1320) mass within er- 
ror bars. The x 2 /d.o.f. of the chiral fit of T2 is larger 
than three (see Tab. IXIljl . where the major contribution 
stems from ensemble A66. Finite volume effects could be 
responsible for the significant deviation of this particular 
value. 



VI. MESONS WITH STRANGE VALENCE 
QUARKS 

In 2-flavor simulations, strange hadrons can be studied 
by including the strange quark just as a valence quark. 
The corresponding quantum field theory is not well de- 
fined, the probability distribution of physical observablcs 
is not anymore strictly non-negative. Nevertheless, since 
the strange quark is heavy compared to the light, dynam- 
ical quarks, observables can be measured and regarded as 
predictions including systematic errors. We stress that 
even though light hadrons are well defined in 2-flavor sim- 
ulations, they also show the systematic error of neglect- 
ing strange sea quarks when the results are compared 



to experiment. From this point of view, the predictive 
power of strange valence hadrons is not significantly be- 
low the one of light hadrons in 2-flavor simulations. The 
strange quark mass parameter is set in each ensemble 
such that the 17(1672) is reproduced (always assuming 
that ruezp = 0.48 fm) (see Sec. MI Bp . 

In contrast to isovector light mesons, C-parity is no 
good quantum number for I = | strange mesons due to 
the non-degeneracy of the light and strange quark mass. 
At unphysically large pion masses, however, C-parity is 
approximately restored. Our interpolators (see Appendix 
[A)) are constructed such that C-parity is a good quantum 
number in the limit of degenerate quark masses. There- 
fore, by monitoring the eigenvectors of the variational 
method, we can learn about the C-parity content of the 
states. 

Since excited states are always more difficult to deal 
with than ground states, this raises the demands on the 
variational method. In some cases it is therefore sugges- 
tive to separate the channels according to C-parity. At 
our largest pion masses, around 600 MeV, one expects C- 
parity to be almost restored. Approaching the physical 
point, C-parity is violated stronger and stronger, and the 
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FIG. 9. (top) (J = |) 0" (X). (bottom) (/ = §) 0+ (if ): 
The S-wave scattering state ttK for zero and minimum non- 
zero relative momentum is indicated for all ensembles using 
crosses. The chiral fits are omitted for clarity. 



corresponding mixing of interpolators is expected to be- 
come increasingly important. To investigate this mixing, 
we include all possible interpolators in the correlation 
matrix, but we also analyze separately the sectors with 
given C-parity. The advantage of the second approach 
is a clearer distinction of the energy levels, where some 
come in the [C = +1] sector, some in the [C = —1] sec- 
tor. In the combined correlation matrix we see both sets, 
but due to the increased noise, fewer levels can be reli- 
ably determined. We discuss this point in the subsequent 
channels. Our results for the dominant C-parity assign- 
ments agree qualitatively with Q. Here we also discuss 
the corresponding mixing, which is accessible due to our 
lighter pion masses. 



Scalars 



FIG. 8. 2 ++ (02) in both representations. 



(I = (K): In the strange (I = channel, 
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interpolator (1) is used for the ground state, which ex- 
trapolates close to the experimental kaon (see Fig. [S]). 
The x 2 /d.o.f. of the chiral fit is larger than four (see 
Tab. lXIlTj) . which indicates that due to the tiny statistical 
errors the systematic errors (e.g. of setting the strange 
quark mass) become visible. For the excited state, we 
use the set (1,2,8,17), its linear extrapolation agrees with 
the experimental if (1460) within error bars. Hence we 
can confirm this state (omitted from the summary table 
of [36|). In this channel we use only h interpolators, 
since the signal of the exotic interpolators is too 
weak, and the corresponding energy levels lie too high. 

(I = ^) + (Ko): The strange scalar channel + is as 
peculiar as its light multiplet partners. The if,* (800) 
(also called k) is a very broad resonance (with a width 
of more than 80% of its mass) and is omitted from the 
summary table of (36j due to its unclear nature. 

Using interpolator (13) alone (not shown), the chiral 
extrapolation almost hits the presumed center of the res- 
onance. To apply the variational method, we use the set 
(10,12,13) and include also (1,4) in the basis at small pion 
masses. We observe that at light pion masses the effec- 
tive masses tend to decrease at large time separations, 
which may be a signal for contributions of a scattering 
state. Like in most cases, we choose a large fit range 
(e.g., 8 timeslices in A66). The results are compatible 
with the if q (800) and the if,* (1430), but also with the 
S-wave scattering state ttK (see Fig. The x 2 /d.o.f. of 
the chiral fit of the ground state is larger than eight (see 
Tab. IXIIip . which is again interpreted as indication for 
systematic errors, probably related to scattering states. 
Here we use only ++ interpolators, the signal of the ex- 
otic H interpolators is too weak. 



B. Vectors 



(I = i) 1 (K*): Considering the strange J p channels 
as mixing of J p+ and J p ~, one can use information 
from the corresponding light J PC channels to speculate 
about the dominating C-parity in the low-lying states 
of the strange J p channel. Based on that analogy, in 
the scalar channels one expects dominance of positive C- 
parity, which is confirmed by our results. In the vector 
channels, however, both C-parities are expected to con- 
tribute to the measurable low-lying states. Looking at 
the experimental states in the corresponding light meson 
channels p(770), tti(1300), p(1450) and p(1570 or 1700), 
one expects that the if* (892) is an (almost) pure 1~~ 
state, while mixing could become important for if *(1410) 
and If* (1680). 

We first discuss sets of purely negative C-parity inter- 
polators. Taking interpolators (1,8,12,17,20), we extract 
a ground state and up to two excitations. The chiral 
extrapolation of the ground state hits the experimental 
if* (892) nicely (see Fig. [TU|) . which is clearly an (almost) 
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FIG. 10. (/ = |) 1" (if*). Results for interpolators re- 
stricted to subsets with [C ~ — ] are shown on top and [C ~ +] 
in the middle. Note that the ground state is missed in the 
[C ~ +] subset. Results with both types are shown at the 
bottom. 
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FIG. 11. 1 + (ifi). Results for the energy levels are shown at 
the bottom. The corresponding eigenvectors for the ground 
state and the first excitation for the lightest pion mass (A66) 
are shown on top. Interpolators (1,2,17) have [C ~ +], (6) 
has [C ~ — ]. Note the dominance of positive (negative) C- 
parity in the ground state (first excitation). Note furthermore 
that there is some mixing in both states, which is allowed by 
the breaking of C-parity towards light pion masses. At our 
largest pion masses, this mixing is suppressed strongly. 



pure [C ~ — ] state. The excitations are a bit high com- 
pared to the experimental if* (1410) and the if*(1680). 

Considering only 1 ^ interpolators, the chiral extrap- 
olation hits the if* (1680). This suggests that mixing is 
important at least for the if* (1680). 

Finally, taking the set (1,8,9,12,16,20,21), both types 
of C-parities are included in the variational method. In 
this analysis, the three lowest states are dominated by 
[C ~ — ] interpolators, where even for the excitations 
the mixing is compatible with zero. A slight mixing is 
observed in ensemble A66, however, the signal is very 
weak, and the corresponding energy levels cannot be ex- 
tracted reliably. One might wonder why we do not see 
a significant contribution of [C ~ -f ] interpolators to at 
least one of the excitations. A possible interpretation 
is that the mixing is indeed weak in this channel at all 
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FIG. 12. (I = |) 2" (if 2 ), representation T 2 : The eigen- 
vectors for the ground state and the first excitation for the 
lightest pion mass (A66) are shown. Interpolators (2,5) have 
[C ~ — ], (6) has [C ~ +]. Note the dominance of posi- 
tive (negative) C-parity in the ground state (first excitation). 
Note furthermore that there is significant mixing in both 
states, which is allowed by the breaking of C-parity towards 
light pion masses. At our largest pion masses, this mixing is 
suppressed. The mixing pattern is similar in representation 
E (not shown). 



simulated pion masses and that there is a further state, 
dominated by [C ~ +], which is not clearly identified in 
the full analysis. The chiral extrapolations of the excita- 
tions come out a bit high compared to the experimental 
if* (1410) and if* (1680), suggesting that simulations at 
smaller pion masses and with higher statistics are neces- 
sary in order to reliably describe the mixing of different 
C-parities and to be able to obtain the if* (1410). 

(I = i) l + (Ki): Looking at the experimental states 
in the corresponding light meson channels ai(1260), 
&i(1235) and ai(1640), mixing is expected already for 
the lowest states if i (1270), if i (1400) and if i (1650). 

Employing pure [C w +] sets of interpolators, the chi- 
ral extrapolation of the ground state ends up between the 
if i (1270) and the if i (1400). The first excitation hits the 
if i (1650) within error bars. From pure [C w — ] inter- 
polators only a ground state can be extracted, the chiral 
extrapolation of which agrees with the if i (1400). 

Allowing for both types of C-parity, three states can 
be extracted when the set of interpolators is optimized in 
each ensemble. The chiral extrapolations are compatible 
with ifi(1270), ifi(1400) and ifi(1650) (see Fig. [IT). 
Since the splitting of if i (1270) and if i (1400) is rather 
small, it is hard to make a statement about its increase 
towards smaller pion masses. (Notice that an increased 
splitting is observed when mixing both charged conju- 
gations for the analogous mesons in the charmed me- 
son sector (Hi)].) This is worsened by the fluctuation 
of the plateau points. However, the eigenvectors indeed 
show stronger mixing approaching the physical point (see 
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FIG. 13. Upper panels: (7 = 2 (7G) in both representations T2 and 75. Chiral fits are suppressed for clarity. 
Lower panels: (7 = |) 2 + in both representations T2 and 77 (reliable signal only in [C ~ +]). 
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Fig. Hip , which is usually accompanied by a more pro- 
nounced splitting. At simulated pion masses, #i(1270) 
and A"i(1650) are dominated by [C « +], #i(1400) by 
[C w — ] interpolators. Our results confirm the existence 
of if 1 (1650) (omitted from the summary table of [36j]), 
which is dominated by positive C-parity in our analysis. 



C. Tensors 



(I = 5) 2~(K2): In the spin 2 channels, investigation 
of the mixing becomes more complicated, since the sig- 
nal is often weak already for the ground state. From 
the light meson states 7^(1670), 7T2(1880) and the (not 
established) p2(1940), one could expect a dominance of 
[C ~ +] interpolators in the ground state. So far, 
#2(1580) is omitted from the summary table of |3(|, the 
lowest established states in this channel are #2(1770) and 
#3(1820). 

Restricting the basis to negative C-parity, we use inter- 
polator (2) as in the corresponding light channel. In both 
T2 and E, the chiral extrapolation is compatible with 
#2(1770) and #2(1820). For positive C-parity, using in- 
terpolator (6) in T 2 arid (8) in E, the chiral extrapolations 
are again compatible with 7^(1770) and #2(1820). 

To take into account both C-parities, the set (2,5,6) 
(resp. (3,4,5,6) in C72) is chosen in T 2 and (2,5,8) in E. 
The two lowest eigenvalues are very close and have to be 
sorted according to the eigenvectors. The eigenvectors of 
T2 are shown in Fig.[T2] We observe that the ground (ex- 
cited) state is dominated by positive (negative) C-parity. 
However, there is significant mixing in both states, which 



appears to be the strongest mixing of all channels con- 
sidered. Strong mixing is also observed in representation 
E. The chiral extrapolations are compatible with the 
experimentally established #2(1770) and #2(1820) (see 
upper panels of Fig.H"3f and do not confirm the #2(1580), 
which is omitted from the summary table of [36j. How- 
ever, increasing mixing towards lighter pion masses could 
still change the slope of the chiral extrapolation. 

(1=1) 2 + (K?j): No experimental state is known in 
the light-quark 2 H channel. In the light 2 ++ channel, 
the lowest states are <Z2(1320), a2(1700) and a2(1950), of 
which the latter two are not established. In the strange 
2 H channel the lowest experimental states are #2 (1430) 
and the (not established) #£(1980). 

The signal of negative C-parity interpolators is weak 
here, thus we restrict our analysis to positive C-parity 
interpolators. Interpolator (2) (Table IXI[) is used in T2 
and interpolator (2) (Table ILX|) in E to extract a ground 
state mass. In both lattice channels, the chiral extrapo- 
lation hits the experimental 7^(1430) nicely (see lower 
panels of Fig. ITB"]) . 

D. Isoscalar light mesons 



1 (0): In principle, correlation functions of isoscalar 
mesons include connected and disconnected diagrams. 
The low lying isoscalar (f> mesons decay mainly into kaons, 
thus one expects that these states are dominated by 
strange quarks (Zweig rule). Since disconnected dia- 
grams are dominated by loops of light sea quarks, it is 
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FIG. 15. Results for the isovector light meson masses. All 
values are obtained by chiral extrapolation linear in the pion 
mass squared. Horizontal lines or boxes represent experimen- 
tally known states, dashed lines indicate poor evidence, ac- 
cording to [3(|. The statistical uncertainty of our results is 
indicated by bands of 1 a, that of the experimental values by 
boxes of 1 a. In case of spin 2 mesons, results for T2 and E are 
shown side by side. Grey symbols denote a poor xVd.o.f. of 
the chiral fits (see Tab. DOTj) . 
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reasonable to assume that the </> mesons are dominated 
by connected (strange) diagrams. We extract <fi meson 
masses evaluating only these connected diagrams, albeit 
with the systematic error of neglecting the disconnected 
diagrams. We use the same set of operators as in the 
light isovector 1 (p) channel (Sec. IV Bp to extract 
three energy levels. 

The ground state mass extrapolates to a value very 
close to the experimental </>(1020) mass (see Fig. [14)l . 
which confirms our choice of the strange quark mass pa- 
rameters. The extrapolation of the excited states ends 
up significantly higher than the experimental 0(1680). 
Since the first excitation p(1450) in the light isovector 
channel is reproduced nicely, one may conclude that the 
neglected disconnected diagrams play a more important 
role for the 0(1680) compared to the 0(1020). The lattice 
irreducible representation T± couples to continuum spins 
1 and 3 (among others). This is why we also indicate the 



FIG. 16. Same as Fig |15l but for strange mesons (left) and 
isoscalars (right). The strange quarks are implemented in the 
partial quenching approximation. The isoscalars additionally 
suffer from neglected disconnected diagrams. In the isoscalar 
sector we show experimental results only where they have 
a dominant ss content [3(J. Grey symbols denote a poor 
X 2 /d.o.f. of the chiral fits (see Tabs IXTTTI and IXlV)) . 



possible spin 3 state 03(1850) in the figure. This state is 
hit by the extrapolations of the first and second excita- 
tion. However, all our interpolators in this channel have 
a naive continuum limit of spin 1. One of these two levels 
may bend down if disconnected diagrams are included. 

2 ++ (f2): As in the meson channel, the experimental 
decay channels of the isoscalar light meson /2 suggest 
dominance of connected diagrams. We use the same in- 
terpolators as in the isovector 2 ++ (02) channel. The 
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results of T2 and E agree (see Fig. IT4"|) . but their chiral 
extrapolations are in better agreement with the /2(1525) 
than with the /2(1430). The latter needs confirmation 
and is not listed in the summary table of (36j . It is un- 
clear if inclusion of the neglected disconnected diagrams 
would yield the /2(1430) or if the ground state of the 
theory is the established /2(1525). 



VII. CONCLUSIONS 

We presented results for the light and strange meson 
spectrum from two dynamical Chirally Improved quarks. 
Seven ensembles with pion masses between 250 and 600 
MeV were analyzed with the variational method in order 
to extract energy levels for ground and excited states. 
In addition to dynamical light quarks we also included 
strange quarks within the partially quenched approxi- 
mation, fitting the strange quark mass by requiring the 
correct £7(1672) mass. 

Figure [T5] shows our chirally extrapolated results for 
the spectrum of light mesons compared to experimental 
values from J3(| . Figure [JS] contains a similar plot for 
strange mesons (left panel) and isoscalars (right panel). 
The results are in general in good agreement with experi- 
ment. For the strange mesons the good agreement for the 
ground states in the kaon, the K* and <p meson channels 
suggest that these observables are well-reproduced in the 
partially quenched approximation and confirm our choice 
of strange quark mass parameter. As discussed in more 
detail in Sections [V] and I VII we do not see any clear in- 
dications of scattering states, which probably show only 
little overlap with the one-particle interpolators used in 
this work. Exceptions are the strange + channel and the 
light isovector + channel at small quark masses, where 
our signal is also consistent with a two-particle scattering 
state. 

The strange meson channels 1~, 1 + and 2~ have been 
investigated with respect to their approximate C-parity. 
In the 1~ channel, the three lowest states seem to be 
dominated by negative C-parity, while positive C-parity 
was shown to contribute to a state in the vicinity of the 
second excitation. The 1 + channel shows some mixing of 
different C-parity towards light pion masses, and the low 
lying spectrum seems to contain states with alternating 
C-parity dominance. The 2~ channel shows strong mix- 
ing towards light pion masses and the ground state (first 
excitation) is dominated by positive (negative) C-parity. 

For our lightest three pion masses finite size effects 
may play a non-negligible role and their influence on our 
results deserves further attention. A study on larger vol- 
umes is in progress and we will investigate this source of 
possible systematic errors in the near future. The larger 
volume will also be used for the spectroscopy of low-lying 
baryon states, where finite volume effects are expected to 
be more pronounced. 
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Appendix A: Tables of interpolators 

In the tables for meson interpolators (Table ITVl to IXII) . 
the two quark fields are labeled by a and b. These are 
placeholders for light (u, d) or strange (s) quarks. The 
indices n, w and di correspond to the smearings narrow, 
wide and derivative, respectively, ji are the spatial Dirac 
matrices, 74 is the Dirac matrix in time direction, 
is the Levi-Civita symbol, Qijk are Clebsch-Gordon co- 
efficients, where all elements are zero except Qui = -^=, 

Ql22 = ~73> ® 2n = ~7e' ® 222 = ~7e and Q233 = 7T 
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13 


ae ; 7,7*756™ +a n j i >yt~t5bd i 




14 


ag^ijfYBbn — a„7i7t75&e i 


+ 


15 


aa^i-ytjabw + a™ 7^74 75^. 




16 


ae i 7i7*756ti' - a™ 7^7*75^ 


+ 


17 


aa 4 75694 


+ 


18 


ae i 7*756s i 


+ 



TABLE IV. Meson interpolators for J p = 0". The first row 
shows the number, the second shows the explicit form of the 
interpolator. In the last column the C parity is given, which 
is only an approximate quantum number in the case of dif- 
fering quark masses. Interpolators with different quark field 
smearings and similar Dirac structure are grouped and these 
groups separated by white space. 



interpolator(s) C 

1 a n -y k b n - 

2 anjkbw + a-wjhbn 

3 a n jkb w — awjkbn + 

4 a w -y k b w - 

5 a n ykjtb n - 

6 a n 7 fe 7t b w + a w 7 fe j t b n 

7 a.n'fkltbw — a w jk7tb n + 

8 a w 'fk'ytbw - 

9 ag k b„ + a n bg k + 

10 7fa fc b 7 i — 

11 ag k b w +a w bg k + 

12 ag k b w — a w bg k — 

1 3 a 9fc 7t b n + a n j t bg k 

14 ag k ^ t b n — a„jtbg k + 

1 5 ag k Jthw + a w -f t bg k 

16 aa k jtb w — a w jtbg h + 

17 ag i -ykbg i - 

18 «s i 7fc7* & a l - 

19 ag k e ijk jj-y 5 b„ + a„e ijk 'yj'y 5 bg k + 

20 ag k e i]k ^j^ 5 b n - a n e ijk 7,75 &e fc 

2 1 a Sfc e 4 jkJjjBbw + a w €ijkJjjBbg k + 

22 ag k e ijk -yjj 5 b w -a w eij k jjjsba k 

23 ag k eijk'jjjt-f5bn + a n eiik'yj"ft'Ysba k 

24 a efc e ijfe 7 j 7 t 7 5 fe n - a„e 4jfe 7j7 t 7 5 fe Sfc + 

25 ag k e i:i kjjjtj5bw + a Ii .eijfc7j7t75&e fc 

26 ag k eijkjjjtjBb w —a w eijkjjjt'YBba k + 



#0+ interpolator(s) C parity TABLE VI. Same as Tab. IIV[ now for J p = 1 . 

I flnfrn + 

3 QjibiD CLw^n 

4 Sm&tu + 

5 ag^ibn + a n jiba i - 

6 ag^/ibn -a„jibg i + 

7 agtfibw +a w jibg i - 

8 ag^ibw -a m r y i bg i + 

9 ag^ijtbn + a„ji-ytbg i - 

10 ag^ijtbn — On 7i 7t + 

II ag^ijtbw +a w ~f i -f t bg i - 

12 ag^ijtbw - a w -fi r y t bg i + 

13 Sa i 6 Si + 



TABLE V. Same as Tab. [TV] now for J p = 
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#1+ 


interpolator(s) 


C 


#2+E 


interpolator(s) 


C 


1 




+ 


1 


Qijkaa k l]b n + Q l jka, n ^ ] ba k 




2 


ttn7fc75&iD + <7™7fc75&n 


+ 


2 


OJ ft, ^ (J k I J 1 L <*J n. f k I J ^ (J k 


+ 


3 


I n. 1 O ^ tu ™uj / n, ] O ly ft 




3 


O ) 7t Oki 4- CJ i 'il-Qjin'y ibpi, 

^6 l 3 ti ' &k 13 w 1 *J re /J^Ofc 




4 




+ 


4 


CjinhCbri, T-i Oil! — C) i <i l- Clin" 1 / i bfi, 

LJft> lJ k 1 J UJ I'jft' I J 


+ 


5 




+ 


5 


Qijkaa k jjjtb n + Q ! jfcan7i7i^e fc 




6 


Tin "y<\bn — a^'y^bn 




(i 


(DijkCiA TiTf^ — (Di ■n-cLn'y'i'ytbf) 

^"C IJ ft' 'J ' " vs 0%jK^ A, 'ii Ij li^Ofc 


+ 


7 


O-d^^bw + &w'YhbQ k 


1 
I 


7 


Qijk<ia k JjJtb w + Qijk<i-w"fj"ftba k 




8 


Tlri, 'Y^&m — (Im'y^bf), 

k J O ^ LU UJ ] O (J fa 




8 


CJi ihdr), 'Yt'Y+Oim O J l' ft 1, i "/ i "/+ Or) . 

^lJK^O k I J JZ^W ^CIJK^W I J \l^O k 




9 




+ 


9 


QijkO-ajba k 


+ 


10 


as fc 7t75&n — a7 l 7t75&s fc 










11 


aa k Jtj5b w 4-a^7i7 5 6 Sfc 


+ 


10 


Qijkaajitbak 




12 


ffle fc 7t75^ -a™7i75bs fc 
















TABLE IX. Same as Tab.|TVl now for J p = 2+£. 




13 




+ 








11 


e-ijkaa^jbn + £ijka n yjba k 











15 


tijkaa k "ijbn — €ijka n 'yjba k 


+ 








16 


(-ijkaa k ^jbw + eijka w -fjba k 


_ 








17 


(-ijkaa k ^jb w — iijka, w ^ijba k 


+ 














#2-T, 


interpolator(s) 


C 


18 


^ijk Ui d k /j /i^ra /j l tU o k 




1 


\£ijk\aa k ljTob n + \£ijk\a n JjTobd k 


+ 


19 


eijkad k jjjtb n - Ujkanjjjtbdk 


+ 


2 


\£ijk\aa k Jjj5b n - \e ijk \a n jj 75 b dk 




20 


£ijkQd k Jj7tbw + £ijkO>w"f3ltbd h 





3 


kijfe |oa t .7j756 w + e ijfe |a ro 7j 75 6e fc 


_l_ 
\ 


21 


Ujk^d k ^j^tb w — eij k a w ^j^ t bd k 


+ 


4 


|eijfe|aa fc 7j75&u. - |eijfe \a w ^ 5 ba k 




22 






5 


|ei 3 fe|as fc 7j7i7s&n + kijfe |a„7j7t7s&S fc 




23 


«n7fc 7*75^ 4-a™7fc7*75&n 




6 


|eijfe|a afc 7j7t756n - |e„fe|a„7j7i7 5 6s fe 


+ 


24 


fln 7fc7t75&™ — a™7fc7*75&n 


+ 


7 


e.j* ae fc 7j7*75&t« + eij/b «™7j7t7s&e fc 




25 


au.7fc7t75&™ 




8 


|eyfe|ae fc 7j'7i756w - |eyfc a«.7j'7t75&e fc 


+ 



26 ae i 7fc7t7s&s l - 

^^^^^^^^^^^^^^^^^^^^^^^^^= TABLE X. Same as Tab. Mi now for J p = 2~T 2 . 



TABLE VII. Same as Tab.HV] now for J p = l+. 



#2-B 


interpolator(s) 


C 


1 


Qijkaa k ~li 7*75 b n + Qijka n 7jlt75ba k 




2 


Qijkaa k ^j r ffysb n - Q i jka„^jy t 'y5ba k 


+ 


3 


Qijkaa k -yj^t'f5bw + Qijka, m ^j^t , y5ba k 




4 


Qijkaa k jjjtj5b w - Qijka w ~fjyt~f5be k 


+ 


5 


Qijkaa j 'y5ba k 


+ 


6 


QijkO,a j 'yt-y5ba k 


+ 


7 


Qijfefla fc 7j75&n + Qijka„jjj5ba k 


+ 


8 


Qijkaa k ^j^b n - Qijfcan7j75&s fc 




9 


Qijkaa k ljihb w + Qijka w ^j^i,ba k 


+ 


10 


Qijkaa k jj't5bw — Qijka-wfjj5ba k 





TABLE VIII. Same as Tab.[lV] now for J p = 2~E. 



#2+ t 9 intcrpolator(s) C 

1 \^ijk\o.d k 7jb n + \^ijk\a n 'yjba k 

2 \tijk\aa k ^jb n — \eijk\a n ^jba k + 

3 \tijk\aa k ^jb w + \e ijk \a w -yjba k 

4 |e l jfc|as fc 7j^ - \^ijk\a w ^jba k + 

5 |eijfe|aa fc 7j7t & n + I a n7j7tfea fc 

6 |eijfc|ffle fc 7j7t & n ■ — \^ijk\a n 'yj'ytba k + 

7 jeijfe|ae fc 7j7t&™ + |a tu 7j7toa fc 

8 |eijfc|ae fe 7j7ifr™ — \ujk\a w ^jj-y t ba k + 



TABLE XL Same as Tab.[TV] now for J p = 2+T 2 . 



light meson 


energy level [MeV] 


X 2 /d.o.f. 


o-+ 


1407(103) 


8.25/5 


0++ 


976(70) 


12.38/5 


0++ 


1689(103) 


6.70/4 




772(13) 


4.65/5 


1~ 


1528(115) 


2.79/5 


1~ 


1733(143) 


2.51/4 


1-+ 


1370(260) 


3.78/5 


1+- 


1347(26) 


8.12/5 


1++ 


1238(33) 


9.62/5 


1++ 


1754(103) 


4.91/5 


2— (T 2 ) 


1849(222) 


3.77/2 


2-{E) 


1965(183) 


2.16/3 


2-+(T a ) 


1745(96) 


5.06/5 


2"+(£) 


1889(139) 


8.96/4 


2++(T 2 ) 


1399(66) 


16.51/5 


2++(£;) 


1379(60) 


6.03/5 


TABLE XII. Energy levels at the physical point and corre- 


sponding xVd-o.f. for the chiral fits of the isovector light 


meson energy levels reported in this work. 


Sources of large 


X 2 /d.o.f. (> 3) 


are discussed in the text. 




strange meson 


energy level [MeV] 


X 2 /d.o.f. 




509(4) 


20.83/5 


tr 


1434(64) 


6.94/5 


0+ 


884(36) 


41.53/5 


0+ 


1323(81) 


6.92/5 


i" 


896(9) 


7.20/5 


i" 


1633(89) 


1.57/3 


i" 


1919(69) 


1.05/3 


i+ 


1339(20) 


3.90/5 


i+ 


1409(17) 


6.76/5 


i+ 


1709(109) 


2.56/5 


2"(T 2 ) 


1750(54) 


5.31/5 


2" (lb) 


1909(52) 


2.21/5 


2-{E) 


1870(75) 


1.51/4 


2~{E) 


1956(71) 


1.17/4 


2+(T 2 ) 


1452(51) 


7.28/5 


2+(E) 


1392(58) 


5.68/5 


TABLE XIII. 


Same as Tab. IXII1 but for strange mesons. 


isoscalar meson energy level [MeV] 


X 2 /d.o.f. 


1— 


994(8) 


6.51/5 




1857(53) 


7.30/5 




1987(40) 


1.41/5 


2++(T 2 ) 


1581(29) 


12.89/5 


2++(E) 


1578(24) 


7.28/5 



TABLE XIV. Same as Tab. IXIIl but for isoscalar mesons. 
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